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1 Introduction 


The recently achieved progresses in experimental sector related to the charm and bottom 
baryons have provided important clues and motivations for the theoretical studies on this 
area. The necessity for a better understanding of the properties of these baryons such as 
their masses, structures and interactions with other particles have increased the theoretical 
interests on them. Their various properties were studied using different methods. For 
instance their masses were studied in Refs. [1-5] (see also the references therein) via various 
methods such as quenched lattice non-relativistic QCD, the QCD sum rule approach within 
the framework of heavy quark effective theory, the constituent quark model, QCD sum rules 
and a theoretical approach based on modeling the color hyperhne interaction. The Refs. [6- 
19] and references therein provide some examples in which their strong and weak decays 
were studied. 

This work provides an analysis of the strong couplings of the heavy Ab(c) and 
baryons to the positive parity nucleon N/ negative parity nucleon N* and heavy B* / D* 
vector meson. Here by N* we mean the excited A^(1535) nucleon with = ^ . Such 
couplings occur in a low energy regime that preclude us from the usage of the perturbative 
approach. The strong coupling constants are the basic parameters to determine the strength 
of the strong interactions among the participated particles. They also provide us a better 
understanding on the structure and nature of the hadrons participated in the interaction. 
To improve our understanding on the perturbative and non-perturbative natures of the 
strong interaction they can also provide valuable insights. Furthermore, these coupling 
constants may be useful for explanation of the observation of various exotic events by 
different collaborations. Beside these, one may resort to these results in order to explain 
the properties of B* and D* mesons in nuclear medium. The nucleon cloud may affect 
properties of these mesons such as their masses and decay constants in nuclear medium 
due to their interactions with nucleons (see for instance the Refs. [20-25]). Therefore, the 
present study is also helpful to identify the properties of these particles in nuclear medium. 

Here, we calculate the strong form factors dehning the strong vertices A^NB*, AhN*B*, 
TibNB*, Tii)N*B*, AcND*, AcN*D*, TicND* and TicAf*D* in the framework of the QCD 
sum rule [26] as one of the powerful and applicable non-perturbative tools to hadron 
physics. By using we then obtain the strong coupling constants among the 

participating particles. This method has been previously applied to investigate some other 
vertices (for instance see Refs. [6,17,27-29] and references therein). 

The paper contains three sections. In next section, we calculate the strong coupling form 
factors in the context of QCD sum rule approach. Section 3 is devoted to the numerical 
analysis of the results and discussion. 

2 The strong coupling form factors 

In this section we calculate the coupling form factors dehning the vertices among the 
hadrons under consideration using the QCD sum rule method. The starting point is to 
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consider the following three-point correlation function: 

n^(g) =^"J d^x I d^y e^d-y {0\r (My) J^^iO) Mx)) |0), (1) 

where T is the time ordering operator and q = p — p' is the transferred momentum. In 
this equation J* denote the interpolating helds of different particles, Ai symbolizes the B* 
or D* meson, B stands for the A^c) or T,h(c) baryons and N shows the nucleon with both 
parities. 

The three-point correlation function can be calculated both in terms of the hadronic 
degrees of freedom and in terms of the QCD degrees of freedom. These two different ways 
of calculations are called as physical and OPE sides, respectively. The results obtained 
from both sides are equated to acquire the QCD sum rules for the coupling form factors. 
For the suppression of the contributions coming from the higher states and continuum a 
double Borel transformation with respect to the variables and are applied to both 
sides of the obtained sum rules. 


2.1 Physical Side 

For the physical side of the calculation one inserts complete sets of appropriate Ai, B and N 
hadronic states, which have the same quantum numbers as the corresponding interpolating 
currents, into the correlation function. Integrals over x and y give 



(0 I Jn I N{p',s')){0 I J% I Ai{q)){N{p',s')Ai{q) \ B{p, s)){B{p, s) \ Jb \ 0) 

(p2 - ml){p'^ - m%){q‘^ - m^) 

(0 I I N*{p',s')){0 I I Ai{q)){N*{p',s')Ai{q) \ B{p, s)){B{p, s) \ Js \ 0) 

(p2 - m|)(p'" - m%.){q^ - m^) 


( 2 ) 


where ■ ■ • stands for the contributions coming from the higher states and continuum and 
the contributions of both positive and negative parity nucleons have been included. The 
matrix elements in this equation are parameterized as 


(0| J^|iV(p',s')) 
(0| J^|iV*(p',s')) 
{Bb{c){p,s) I I 0) 

(0 I rM I M{q)) 
{N{p',s')Ai{q)\B{p,s)) 

{N\p\s')M{q)\B{p,s)) 


XnUn{p',s'), 
XN*"fbUN* (P^ s'), 
(p, s), 


€’'Un{p',s' 


9ilv 


la. 


va a 

Q 92 


e''U]y*{p',s')75 


9ilu - 


ms + niN 


ub{pa), 


ms + rriN* 


q 92 


ub{p,s), (3) 


where Xm{n*) and Ag are the residues of the related baryons, U]sf(^N*) and are the spinors 
for the nucleon, Afe(Ac) and Sb(Sc) baryons; and /m represents the leptonic decay constant 
of Here gi and p 2 are strong coupling form factors related to the couplings of 
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the B baryon and M. meson to the positive parity nucleon A^; and and g 2 are those 
related to the strong vertices of B baryon and M. meson with the negative parity nucleon 
N*. Application of the double Borel transformation with respect to the initial and hnal 
momenta squared yields 


^ p, mg 

Bn y{q) = Xb/mG ^ e [$17^ + $2 + $3 + 


(4) 


where 

$1 = 


rriM 


„*2 


{niB + - q^) L 

-)- e Atv* 


\N{gi + g 2 ){mB + rriN* 


m% + uiNfUB + 


$2 — 




ms + niN*) + glims - ttin*) ) ( + mM*mB - q^ 

1 

. \ f o-yn 

M 

2 2 


1 r / \ 


$3 — 


+ XN*imB - mN*){glimB + mN*y - glmj^ 

2mM 


{ruM + mN){mB + - q^ 

™‘'n / 

eX^X^yms + mN)yglimB + uin*) - glm^* 
msmM 


eXX^ Xn (ms + ) gi (ms + niN) + g2mN 


0*2 


- eXX^XNigi + g2)imB + tun*) 


imB + mN*)iml^-q'^ 

-j- cXP^ Atv* ( glims A glims — 


(5) 


with and M'^ being the Borel mass parameters. 

2.2 OPE Side 

For the OPE side of the calculation, the basic ingredients are the explicit expressions of the 
interpolating currents in terms of the quark helds, which are taken as 


"^^blc] (^) 

= eabcU^'^ ix)C-i^(fix)ih[c]yix), 


(^) 

= eafec(^^“'^(a:)C'7^d^(a:)j757^(6[c])'’(a;), 


XfNiy) 

= iy)C-ifiuyy)^-f^-iyd\y), 


<Ab*[d*](0) 

= h(0)7^6[c](0), 

(6) 


with C being the charge conjugation operator. By replacing these interpolating currents in 
Eq. (1) and doing contractions of all quark pairs via Wick’s theorem, we get 

[ d*ye~^P"^e^P''yeabceiji 
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X Ur.ll,S‘J(y - xY„CSf(y - x)CyeSt{y)y„St(,^{-x) 

- 757/jS5'’(!/-a:)75CS“’'(!/-2:)C7f>S“(!/)7„Sj"|(-2:)| , (7) 

for Ai,N^*'>B* and AcN^*'>D* vertices and 

J d*ye~'^P"^e^P''yeabceije 

X ^A5AySa\y - x)'y^CSt"^{y - x)C-f0S^^{y)'y^Sl^(y^{-x)-f^'y5 

- 75lyS^J{y-x)-fuCS^^^{y-x)Cj^S^J;^{y)-ff,Sl^(^]{-x)-f^j5^ , ( 8 ) 

for TibN^*W* and T,cN^*^D* vertices. In these equations, <S'^'[c](^) ^u[d](^) correspond 

to the heavy and light quark propagators, respectively. Using the heavy and light quark 
propagators in coordinate space and after lengthy calculations (for details see Refs. [28, 29]), 
we obtain 


n 


OPE 


(q) 


4p, + 4i, 

other structures, 


(9) 


where the 11*(q'^) functions contain contributions coming from both the perturbative and 
non-perturbative parts and are given as 


n 


OPE 


{q^) = ds ds' 


pert/ / 2\ I non—pert / / 9 

X (s,s,g^) + Pi ^ (s,s',g" 


(s — p^){s' — p' ) 


( 10 ) 


The spectral densities pi{s,s',q‘^) appearing in this equation are obtained from the imag¬ 
inary parts of the Ilj functions as pi{s,s',q‘^) = -Im[Ili\. Here, as examples, only the 
results of the spectral densities corresponding to the Dirac structure 7^ for A^NB* vertex 
are presented, which are 


ffr\s,s\q^) = 


nibmuS 

647r^Q 


■0 


Li(s,s',g^) 


f*l—X 


+ 


dx 


dy- 


mlx‘^{X' + 2x) 


647r4T’3 

X {X + y) + q^xy \3y{X’ - + 4x) - 2X'\3x + X') + 2y‘^{l3x^ - 14x + 2) 

— 2gX xyX' (^s(2 + ISa;^ - 18a;) - s\X' - 3)) + (^2sx(l - 10a; -f Ihx^) 

— 4sa;^T’'^ — s'(l — 21a; -|- 41a;^ — 15a;^) j -|- s'y^{l — 17a; -|- 30a;^) — Asx^X'^ 

+ mlx(^mu{3 — 5x + 2x^ — 2xy) + 3md{X' + x)x'^ + X"^ s'^y{5y — 8xX' 

— 3Axy — 6y^ + 15x'^y + 30xy^) + 3s‘^x‘^{l — Ay — 6x + 3x^ + IQxy) + 2ss'x 

X {by — 9y^ — AxX' + Ibx^y — 26xy + 30xy‘^) + 2mlx q^ (^3xX' — 3y + 16xy 

+ 8x^y — Ay'^ + 16a;p^ j -|- X (^s'{3y — 3xX' — 16y + 8xy — Ay"^ + IQxy'^) 


A 



+ 2sx{l — ?>y — 'bx + Ax^ + Sxy)^ + mj, SnidX {^sxX{X' — 1) + y{s'u — q^x) 

X (3y — 1) j + mu {^sxX (6 — 9a; + 3a;^ — 3xy) + y {q^x{Ax — 3x^ + 3xy + y — 1) 


+ Ss'XX'" -3s'yX{X' + 2) 



■0 


L2{.s,s\q ) 


and 


non—pert / / 2\ 

Pi (s,s,g) = 


{uu) 

pl—X 

dx dy 

Jo Jo 


s'{mu - 2mb) - q^md 0 Li{s, s', q^) 
1 




x + X')X — mb{x + XX') — muX 


+ (uu) imu{3xy + 3xX' — y) — mb{x + X') — 2mdX] 0 ^ 2 ( 5 , s', (f) 


G\ 1 

- («.—)■ 


71 ' 11527r2Q4 


9mbQ'^{md - 2m„) + 


+ m^ 


+ 3s'" ( - 2g^mfe(mb - rrin) + 0Li(s,s',g^) 

1 


el—x 


3A’''(2x + A") + y^(l3 + x(39A’' - 


+ y{2x + A") + X'{\\X' - 1) + 6?/3(2x + A") | 0 | L 2(s, s', g') 
1 


1927r2Q 


ml{dd){Qmb + dm^) + ml{uu){7mu — 3md — 18mb) 


X 0 


Li(s,s',g ) 


( 11 ) 


( 12 ) 


where 


X = X + y — 1, 

A" = a:-l, 

Q = ml- q^, 

Li{s,s',q^) = s', 

L 2 {s,s',(f‘) = —m^a; + sx — sa;^ + s'g + g^a;g — SX7/— s'a;g — s'l/^. (13) 


The 0[...] in these equations is the unit-step function. As already stated, the match of the 
results obtained from physical and OPE sides of the correlation function gives the QCD 
sum rules for the strong coupling form factors. As examples, for the form factors related 
to the AbNB* and AbN*B* vertices, we get 


Pi(g^) = 


niiv (m|. - g^) I 2 

e^eM' -—-— {mjd 

Ajy tt 


m\^{Yiy, — 2112) — 2VmArmAr*n4 — 2m\^{mN*B.2 + 114) 
m\,^ ^2g^n4 -|- {mJjq — mNmN*){mM*'B^ + 2114) j 

m\ (m^rnN *(2112 - Hs) -h m^*n3 + m^(n3 - 2112) -h 2mAr*n4 - 2111 
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+ {rriAi, — rn]sf*){mA^ + rn]\f*) mAi,m]\fm]\f*VIls — 2mNrn]\f*VIli + + 2q^mA^Il4^ 

+ ml^(2Ili - {itinV + m%,)Il3] |, 

. 'i r 

- Ila + mTvml Ilg + mA,m%,VIl3 


) = e’M^CM'- 


^ (m|. — q'^)(mAr + ttiaJ 


XnH 


+ mAr»m^^(2mAr»n3 — mArlls — 2114) + 2 mA^ {mN{rn\i, — q^) + mN*q^'^^i — 2m^»Vn4 

+ - V) (^mAt(2mAtn2 - mA^IIs + mAr.IIs) + 2(mA(, - mAr*)n4j 

— m\^ {mN*{mNrriN*^z — ArriN^A + 4mAr*n4) + 2Vnij 


!I!iv (m|. — I , 14 ^ 

= e^cM' —^—<; (mAt - mAr)(mAb + mAf) m* IIs + mA.mArmAT.Vns 




+ 2mArmAr»Vn4 + 2mAtg 114 + ”^A6(2ni — {uinV + m^.)n3 


+ rriB, 


mAifnM*'mA^VI\.‘i + 2 mA^(mArn 2 — 114) — 2mArmAr*Vn4 


+ - 2112) - 2mA, (^mTv.V - + m\^ (^m^*(2n2 - 113) - m^n3 

N 

+ mAT (2114 + mAr» 113 — 2mAr» 112) + 2ni 


:ilNim%*—q^){mA,+mN*) 

= P m2 pm' - 

XN*n 


"^A, + + mA,mAr* jn3 + m\^mN 

X ^(V — mAr)n3 — 2114j — 2m%VIl4 + 2mA, + g^Vj 114 — m|* (mA, — V) 


where 


X I mA, [2mAjl2 - (mA, + mAr)n3 ) + 2(mA, + mAr)n4 
— 4mArn4 + 4mAr*n4^ — 2Vni 


H = 


<6 


mTV ( mArm7V*n3 


(14) 


V = 


2/B*AA,mB*mA,(mA, - V)(m7v + mjv*) (m%, + m\^ - mA^mA, - m% 

mN*{mN + mA,) - m%, 

{niN - niN*). 


(15) 


3 Numerical results 

To numerically analyze the sum rules for the strong coupling form factors and to hud their 
behavior with respect to = —g^ we need some inputs as presented in table 1. Besides, we 
also need to determine the working regions corresponding to the four auxiliary parameters, 
M^, M'^, So and Sg. The and emerge from the double Borel transformation and Sq 
and Sg originate from continuum subtraction. These are auxiliary parameters, therefore, we 
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Parameters 

Values 

{uu){l GeV) = {dd){l GeV) 

-(0.24 ±0.01)3 Qgys [ 3 Q] 

(^) 

(0.012 ±0.004) GeV^ [31] 

1 

0 

(0.8 ±0.2) GeV^ [31] 

mb 

(4.18 ±0.03) GeV[32] 

rric 

(1.275 ±0.025) GeV[32] 

rrid 

4.8l[]i MeV[32] 

rriu 

2.3+°;^ MeV [32] 

ms- 

(5325.2 ± 0.4) MeV [32] 

mo* 

(2006.96 ± 0.10) MeV [32] 

niN 

(938.272046 ± 0.000021) MeV [32] 

TTljv* 

1525 TO 1535 MeV [32] 

rrihb 

(5619.5 ± 0.4) MeV [32] 


(2286.46 ± 0.14) MeV [32] 


(5811.3 ± 1.9) MeV [32] 


(2452.9 ± 0.4) MeV [32] 

fs* 

(210.3l?i) AleV [33] 

/d* 

(241.91}^;}) MeV [33] 


0.0011 ±0.0005 GeV^ [34] 

Xn* 

0.019 ± 0.0006 GeV3 [35] 

Aaj 

(3.85 ±0.56)10-2 GeV3 [36] 

Asi, 

(0.062 ±0.018) GeV^ [37] 

Aa. 

(3.34 ±0.47)10-2 GeV^ [36] 

Aec 

(0.045 ±0.015) GeV^ [37] 


Table 1: Input parameters used in the calculations. 


need a region of them through which the strong coupling form factors have weak dependency 
on these parameters. The continuum thresholds are in relation with the hrst excited states 
in the initial and hnal channels. To determine them the energy that characterizes the 
beginning of the continuum is considered. Table 2 presents intervals of the continuum 
thresholds used in the calculations. To determine the working regions for the Borel mass 
parameters, we need to take into account the criteria that the contributions of the higher 
states and continuum are sufficiently suppressed and the contributions of the operators with 
higher dimensions are small. The intervals obtained based on these considerations are also 
given in table 2 . 

The determination of the working regions of auxiliary parameters is followed by the 
usage of them together with the other input parameters to obtain the variation of the 
coupling form factors as a function of Q^. For this purpose, the following fit function is 
applied 


9bnm(Q^) ^ Cl + C2exp 



(16) 


where Ci, C 2 and C 3 for different vertices are given in tables 3-6. This £t function is used 
to attain the coupling constants at = —nij^ for all structures. The results for coupling 
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Vertex 

SoiGeV^) 

s'oiGeV^) 

M\GeV^) 

M'\GeV^) 

AbN^*^B* 

32.71 < So < 35.04 

1.04 < s[, < 1.99 

10 < m2 < 20 

1 < M'" < 3 

SfoVWR* 

34.91 < So < 37.40 

1.04 < s[) < 1.99 

10 < m2 < 20 

1 < M'" < 3 


5.71 < So < 6.72 

1.04 < s[) < 1.99 

2 < m2 < 6 

1 < M'" < 3 

ScA^WD* 

6.51 < So < 7.62 

1.04 < s[) < 1.99 

2 < m2 < 6 

1 < M'" < 3 


Table 2: Working intervals for auxiliary parameters. 

constants are presented in table 7. The presented errors in the results arise due to the 
uncertainties of the input parameters as well as uncertainties coming from the determination 
of the working regions of the auxiliary parameters. From this table we see that the maximum 
value belongs to the coupling constant associated to the vertex AhN*B* and the minimum 
value corresponds to the coupling gi related to the AcND* vertex. 



Cl 

C2 

C3(GeV2) 

9im 

-2.44 ±0.68 

-0.34 ±0.10 

-17.88 ±5.18 

92{Q^) 

22.92 ±6.64 

3.87 ± 1.12 

16.85 ±4.89 

9l{Q^) 

-6.21 ± 1.73 

-26.76 ±8.01 

-193.72 ±56.17 

92m 

88.27 ±25.60 

9.65 ±2.70 

24.32 ± 7.05 


Table 3: Parameters appearing in the £t function of the coupling form factor related to the 
AhN^*'>B* vertex. 

In conclusion, we calculated the strong coupling constants related to the vertices AbNB*, 
AbN*B*, T.bNB\ T.bN*B\ A^ND*, AcN*D*, E^ND* and in the framework QCD 

sum rules. Our results may be checked via other non-perturbative approaches. The pre¬ 
sented results can be helpful to explain different exotic events observed via different exper¬ 
iments. These results may also be useful in the analysis of the results of heavy ion collision 
experiments as well as in exact determinations of the modihcations in the masses, decay 
constants and other parameters of the B* and D* mesons in nuclear medium. 
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Cl 

C2 

C3(GeV2) 

9im 

297.08 ±89.12 

-282.66 ±81.97 

-225.11 ±67.53 

92 m 

-18.12 ±5.07 

-3.82 ± 1.14 

14.06 ±4.08 

9l{Q^) 

87.60 ± 25.40 

-82.34 ±23.87 

-24.88 ±7.21 

9im 

31.80 ±9.22 

0.90 ±0.26 

-6.70 ± 1.94 


Table 4: Parameters appearing in the fit function of the coupling form factor related to the 
vertex. 



Cl 

C2 

C3(GeV2) 

9im 

1.28 ±0.36 

0.92 ±0.27 

-155.98 ±45.24 

92 m 

3.88 ± 1.13 

1.27 ±0.38 

3.60 ± 1.04 

9im 

3.01 ±0.87 

(17.97 ±5.21)10“^ 

-2.77 ±0.80 

92{Q^) 

11.52 ±3.23 

-2.38 ±0.71 

2.26 ±0.66 


Table 5: Parameters appearing in the £t function of the coupling form factor related to the 
vertex. 
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